Three-Body Recombination in One Dimension 
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We study the three-body problem in one dimension for both zero and finite range interactions 
using the adiabatic hyperspherical approach. Particular emphasis is placed on the threshold laws 
for recombination, which are derived for all combinations of the parity and exchange symmetries. 
For bosons, we provide a numerical demonstration of several universal features that appear in the 
three-body system, and discuss how certain universal features in three dimensions are different in 
one dimension. We show that the probability for inelastic processes vanishes as the range of the 
pair-wise interaction is taken to zero and demonstrate numerically that the recombination threshold 
law manifests itself for large scattering length. 



I. INTRODUCTION 

One-dimensional (1-D) few-body and many-body sys- 
tems have been the subject of intense theoretical study 
for many years [l|, 0, H, Q, H, 0| . This is largely because 
certain one-dimensional problems admit exact solutions 
using the Bethe ansatz. These theoretical studies are 
gaining increasing attention due to the experimental re- 
alization of effective 1-D geometries in tightly confined 
cylindrical trap geometries @, H, H, 0, EL G2 ■ Atoms 
in such traps are essentially free to propagate in one co- 
ordinate while being restricted to the lowest cylindrical 
mode in the transverse radial coordinate. Olshanii has 
shown [l3[ how the 1-D scattering length (which we call 
a) is related to the three-dimensional s-wave scattering 
length a?,D and to the oscillator length a± in the confined 
radial direction. This identification allows a connection 
with the intensively-studied 1-D zero range model with 
two-body interactions 



V(x) = gS(x). 



(1) 



The coupling constant g is renormalized to account for 
virtual transitions to excited radial modes which appear 
as closed channels in the unrestricted coordinate x. Early 
analytic calculations by McGuire [l[ using the interaction 
Eq. (QJ gave vanishing probability for all inelastic events 
such as collision induced break-up (AA + A — > A + A + A) 
and three-body recombination (A + A + A — > AA + A) . 
Section llV Bl discusses how finite-range interactions break 
the integrability of zero-range models, and also how the 
probability for recombination and break-up behave with 
respect to the range of the interaction. In Section IVI1 
we comment briefly on the relevance of our 1-D model to 
physical systems in actual atomic waveguides. 

In three dimensions (3-D), universal features occur 
when the scattering length is the largest length scale 



(see [3] and references therein). Adiabatic hyperspheri- 
cal studies in three dimensions have provided a great deal 
of insight into such universal behavior [l^,[3[I3,[3- For 
example, the Efimov effect in 3-D appears as a hyperra- 
dial potential curve in the region r$ <C R <C |a| that is 
attractive and varies as R~ 2 . Since it has a supercritical 
coefficient, this potential has (as a-^D oo) an infinity of 
long-range bound states spaced geometrically by a uni- 
versal constant [lj| . One of our goals here is to determine 
what kind of universal behavior, if any, appears in one 
dimension. We discuss these issues in Section IVl 

Finally, a major portion of this work deals with the 
threshold laws for recombination in 1-D. We outline how 
these laws can be extracted from the asymptotic form of 
the adiabatic potential curves with a generalized Wigner 
analysis [2(J El, [22j , and show that — as with any such 
analysis — they are independent of the short-range prop- 
erties of the interactions. This work closely parallels the 
3-D analysis we carried out in Ref. [l^j. Threshold laws 
are found for all combinations of the parity and exchange 
symmetries, including the cases where only two of the 
three particles are identical. In order to demonstrate the 
threshold behavior, we present numerical calculations for 
bosons. 



II. HYPERSPHERICAL COORDINATES 

Since hyperspherical coordinates and the adiabatic hy- 
perspherical representation play a central role in this 
paper, but may be unfamiliar to some readers, we will 
briefly discuss the important points. For more de- 
tailed information, we refer the interested reader to 
Refs. dEEllll- 

We separate the center of mass motion from the rela- 
tive motion using Jacobi coordinates, 
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X12 — x\ — x 2 and x 12 ,3 



miXi + m 2 x 2 
mi + m 2 



a*. (2) 



The positions Xi locate each particle relative to some 
laboratory-fixed origin, and m, are their masses. The 



2 



Jacobi coordinates X12 and X12.3 constitute a Cartesian 
coordinate system. Transforming these to polar coordi- 
nates (R, (f>), 



HR 2 = [i\2x\ 2 + A f i2, 3^12,3 an d tan</> = 



Ml2,3 £12,3 



M12 Z12 
(3) 

gives the hypersphcrical coordinate system. The reduced 
masses are 

mim 2 , (mi+TO 2 )m 3 

M12 = ; and yUi2, 3 = ; ; , (4) 

mi + m 2 mi + m 2 + TO3 

and we choose the three-body reduced mass /1 to be 

- v /Ml2Ml2,3- 

With the above definitions, we can find the interparti- 
cle distances to be 



\xi - Xj\ = dijR\sin(<f) - 4>ij)\, 
where we have introduced the constants: 



(5) 




(0) 



The coalescence points — where the interparticle dis- 
tances are zero — are defined as 



tan < 



>2.; 



mi //ii2,3 



m-12 V M12 



tan 03i = 



m 2 //ii2,3 



TO12 V M12 

If all three particles have the same mass, then dij 



(7) 



2/V3 and </>23 = — </>3i = 7r/6. If there are only two 
identical particles, then it is convenient to label them 1 
and 2. In this case, d 2 % = ^31 and 023 = — 031- 

Since the hyperradius R is the only length scale in 
the system, giving the overall size of the system, it is 
natural to treat it as a "slow" coordinate for an adiabatic 
representation [24[ . The Hamiltonian for the system can 
be written as 



where 



(R, 



1 d 2 
2(iR 2 d0 



+ V(R,4>) (9) 



and V includes all of the interactions. Atomic units have 
been used here and will be used throughout this work. 



The adiabatic representation is then defined by the equa- 
tion 



H ad (R,4>)$ u (R] 



U V {R)<S>„{R; 



(10) 



The eigenfunctions & u (R;(f)) are the channel functions, 
and the eigenvalues U y {R) the potential curves corre- 
sponding to each channel. 

In the limit R — > 00, the potentials approach either the 
energies of the bound diatomic molecules for the recombi- 
nation (two-body) channels, or zero energy for the three- 
body continuum channels. In this limit, the channels are 
uncoupled, although there must be coupling at smaller R 
for any inelastic transition such as recombination to oc- 
cur. For the purposes of determining the threshold laws, 
however, we need only know that the channels in this 
representation become uncoupled asymptotically. 



III. THRESHOLD BEHAVIOR 

At ultracold temperatures, the dominant character of 
a given scattering process is controlled by its threshold 
behavior, and the adiabatic hypersphcrical picture read- 
ily yields this behavior [22j ■ We show this by first solving 
the adiabatic equation Eq. (fit))) while taking into account 
the appropriate exchange symmetry and parity. In the 
limit R — > co, all of the adiabatic potentials take the 
general form U K k 2 / (2{iR 2 ), yielding the hyperradial 
equation 



RdR\ dR 



— - h 2 



F(R) = 



(11) 



where k is related to the total energy E (E—0 at the 
three-body breakup threshold) by k 2 = 2fiE. The final 
momentum in the two-body channel is related to the total 
energy by k 2 = 2fi{E + B2), where B 2 is the two-body 
binding energy. This equation is simply Bcssel's equation 
with general solution 



F K (kR) = AJ K (kR) + BY K (kR). 



(12) 



The coefficients A and B are determined by the usual 
procedure of matching to short-range solutions and are 
related to the S'-matrix. It is more convenient for the 
present discussion, though, to consider each S'-matrix el- 
ement for recombination at a fixed total energy E in the 
form 



S fi (E) ~ (F2B®2b\S\F 3 b&3B 



(13) 



where 2B labels the final atom-dimer channel and 3B 
labels the initial three-body continuum state. To deter- 
mine the threshold behavior, we use the small argument 
form of the Bessel functions. For recombination, the fi- 
nal momentum in the two-body channel is non-zero and 
slowly varying at the three-body threshold, but the ini- 
tial momentum k vanishes there. Hence, it is the ini- 
tial channel that determines the energy dependence of 



3 



the transition probability, and it is the lowest three-body 
continuum channel, with k = /€ m ; n , that will dominate at 
threshold. As a result, in the limit k — > 0, the recombi- 
nation probability must scale like 



\S 



fi\ 



(14) 



With this scaling, we can, of course, determine the scaling 
behavior of the recombination rate. This connection will 
be discussed in Sec. IIV Al Moreover, when the scatter- 
ing length is the largest length scale, simple dimensional 
arguments (namely that a probability must be unitless), 
imply that 



|S/*| 2 oc(fca) : 



(15) 



A. Three identical bosons with 5-function 
interactions 

In this section, we will use ^-function pair potentials, 
Eq. ([l}, to find K m ; n for three identical bosons. Even 
though the recombination rate for such interactions ac- 
tually vanishes [if, the analytic solutions of the adia- 
batic equation possible with these interactions [25|, [26[ 
nevertheless give K m in for general interactions. With 
5-function interactions, we can treat only interacting 
bosons. The case of interacting fermions will be con- 
sidered in Sec. IIIIBI 

For three identical particles, the coalescence points 



(e.g. Xi 



0) form radial lines that are equally spaced 



by 7r/3 in the two-dimensional space spanned by the Ja- 
cobi coordinates defined in Eq. (|2|). These lines divide 
the coordinate space into six regions, each of which cor- 
responds to a unique ordering of the three particles along 
the real line. Symmetry thus permits us to solve the adi- 
abatic equation (| 1 0[) in just the region < <p < ir/6 with 
appropriate boundary conditions (see App. [A"|) . 

The (5-function coupling constant g is related to the T 
D two-body scattering length by a = —1/ ((J.2B9) which is 
defined from the 1-D effective range expansion in the even 
parity "partial wave", fctan(<5) = ± + ?fk 2 + ... [H,[13|. 
We require g < so that the potential supports a two- 
body bound state, and write the general solution for this 
channel as 

®(R;<j>) = A sinh q<p + B cosh qcf> (16) 
where q is related to the potential energy through: 



U(R) = - 



(17) 



For even parity, we impose boundary conditions $'(0) = 
and 

lim $'(</. + e) - $'(0 O - e) = 2^gaR^ ), (18) 

where a ~ \jdyi = 3 1 / 4 2~ 1 / 2 and 4>o = n/6. We require 
the additional symmetry condition that <f>'(7r/6 + e) = 



— $'(7r/6 — e). These conditions lead to a transcendental 
equation for q: 



gtanh— = —/igaR. 
6 



(19) 



A similar analysis for the continuum solutions begins 
with the general solution 

$(i?; <j>) = A sin(K0) + B cos(k0) (20) 

where n is related to the potential energy through 



U(R) 



2fiR 2 



(21) 



and is the solution to the following transcendental equa- 
tion: 



K7T 

k tan — = [igaR. 
6 



(22) 



As R — > 00, the allowed solutions are k = 3, 9, 15.... Note 
that Eq. (flT))) is, of course, the analytic continuation of 
Eq. ([22]) with k —> iq. 

For odd parity solutions, the only difference is that 
the boundary condition at <fi = changes to $(0) = 0, 
immediately giving B — and leading to 



and 



q coth — = —figaR 
6 



k cot — = agaR. 
6 



(23) 



(24) 



From Eq. (|24|) . the allowed values k for odd parity bosons 
asymptotically are k = 6, 12, 18 . . . Note that the even 
and odd parity potential curves for the two-body channel 
become degenerate in the limit R — > 00 as one would 
expect. 



B. Three identical particles with square- well 
interactions 

In order to determine whether the 5-function results of 
Sec. lIII Al are general, we now consider square- well pair- 
wise interactions: 



V(xi_ 



-V , if0<|ztf|<L; 
0, otherwise. 



(25) 



Like the (5-function interactions, the adiabatic equation 
remains analytically solvable for square- well interactions. 
In addition, square- well interactions are a good model for 
short — but non-zero — range interactions. At the end of 
this section, we will point out that the results are, in fact, 
general for all short-range interactions. Unlike the 8- 
function interactions, though, indistinguishable fermions 
can interact via square-well interactions. Consequently, 
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We thus focus on solutions that are either completely 
symmetric or completely antisymmetric — assuming as 
in Sec. MI Al that the particles are spin polarized such that 
the spatial wave function carries all of the permutation 
symmetry. As in Sec. 1111 Al symmetry permits us to to 
solve the adiabatic equation (flO|) in a wedge of width 
■|. To simplify the imposition of boundary conditions, 
though, we choose the interval § < < ^f- so that the 
edge (j)b of the square- well centered at <f> = ^ is simply 



COS (f>b 



L 
R' 



(26) 



This condition is invalid at small R where the two-body 
interactions overlap, but is sufficient for determining the 
allowed K in the region R 3> L. Again using U{R) = 
K 2 /(2fiR 2 ) (21]) and defining /3 2 = 2^R 2 Vb, the adiabatic 
equation becomes: 



d 2 



d 2 



< 6< 



< 6 < 



2tt 



(27) 



We write the continuum solutions as 



A sin 



B cos 
C sin[K(^ 



V^ 2 

_ 27T ^ 

3 / 



|--Dcos [k((, 



27T ' 

3 / 



6h < 



< , 



■2k_ 

3 



A similar expression may be written for the two-body 
channels, but we want to focus on the recombination 
threshold behavior and thus need only the asymptotic 
behavior of the three-body channels. 

For brevity, we sketch the derivation only for even- 
parity bosons and summarize the results for all other 
symmetries in App. [Bj For this symmetry, we impose 
boundary conditions: 



Since j3 oc R, this equation has the same form as Eq. (|22p . 
with the same consequences. In particular, the allowed 
k are 3, 9, 15.... 

This analysis can be easily extended to find all allowed 
values of k for all symmetries (see App. [Bj. The thresh- 
old law, however, depends only on K m i n . Collecting this 
quantity for each symmetry, we find for bosons that 

J 3 Even Parity 
Kmin = 1 6 Odd Parity ' 



in agreement with the results from Sec. IIII Al for <5- 
function interactions. For fermions, we find 



6 Even Parity 
3 Odd Parity 



(32) 



These fermion K m ; n are actually the same as one would 
find by symmetrizing the free particle solutions [see 
Eq. (|A6|) ]. The interacting boson K m i n above, though, are 
not the same as the symmetrized free particle solutions 
from Eq. (|A5|) . From Eq. (fl~4"|) , we sec that recombination 
at threshold will be dominated by the even parity symme- 
try for bosons and by odd parity for fermions. Moreover, 
recombination of bosons and fermions (28l | in 1-D share 
the same threshold law. 

The above analysis, in fact, generalizes to short-range 
potentials of any form. That is, if we take the log- 
derivative b [see Eq. (|C3[) ] from the two-body equation 
outside the range of the interaction, then the matching 
condition is 



— ) bR = 



M0b - —j 



(33) 



for even-parity bosons. In the limit R — > oo, b may be 
regarded as a constant so that this equation reduces to 
the same form as Eq. (|22|) or Eq. (|3 1 [) and carries the 
same consequences. 



2ti\ 



(28) 



which gives A = C = 0. 
of the wave function at cj. 
condition: 



Matching the log-derivatives 
= fib yields the quantization 



^n 2 + f3 2 tan ^ K ^+/3 2 (fi b - - ) 



Ktan 



U 2n \ 



(29) 



In the limit R — » oo, fit — » -|, y 1 n 2 + (3 2 — * (3, and 



V« a +/?*(&--)-► VWoh) L = X. (30) 



(31) 



So, in this limit, Eq. (|29j) becomes 



/3tan A = — Attan ■ 



C. Two identical particles with 5-function 
interactions 

For this case, all combinations of parity and exchange 
are considered assuming the spins of each particle are 
fixed by spin polarization. Following the discussion in 
SecHU we label the two indistinguishable particles 1 and 
2. The interaction then takes the form 



V = gs5{x 12 ) + g D [S(x 31 ) + 6(x 23 )] 



(34) 



where gs denotes the same-particle coupling, go denotes 
the different-particle coupling, and are the interpar- 
ticle distances. The constant a in Eq. (fT5|) is equal to 
1/dij for each pair of particles ij. We shall use as for 
particles 1 and 2 and qd otherwise. 

By symmetry, we need only solve the adiabatic equa- 
tion in the range < < 7r/2. The (5-function boundary 
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condition Eq. {18]) is imposed at 0o — 023 for the dis- 
tinguishable pair [see Eq. 0] and at O = tt/2 for the 
indistinguishable pair. The general continuum solution 
is now conveniently written as 

\ Acos (k4>) + B sin (k4>), < < 023; 

= |CcOs[K(0-f)]+-Dsin[K(0-f)] ! 023 < < f ■ 

.(35) 

We are now prepared to determine the asymptotically 
allowed values of k for specific symmetries. 



For brevity, we outline the derivation only for even- 
parity bosons and summarize the results for all cases in 
App.El As before, even-parity bosons require the bound- 
ary conditions 

$'(0) = and *'(- + e) = -*'(-- e). (36) 

These conditions immediately give B — and yield the 
quantization condition: 



-K H — u 2 g 2 asaDR 2 

K 



tan 



6 23 - ^) 



— [igasR tan K023 tan 
I 



»23 



7T, 

2" 



Attan K023 = M3(as + 2a£>)R. (37) 



In the limit i? — » 00, the i? 2 
equation reduces to 

- tan k (023 



term dominates, and this 



2' 



0, 



for n 



(38) 



1,2,3... 



yielding solutions k = n7r/|023 — §| 
Careful inspection of Eq. (|37|) leads to the additional so- 
lution 



k tan K023 - *■ 00 • 



(39) 



This equation has solutions k = (n + l/2)7r/|023| for n = 
0,1,2.... 

Similar analyses can be carried out for all combina- 
tions of parity and exchange (see App. lBj) . The resulting 
Kmia f° r each symmetry are summarized below. As dis- 
cussed above, n m in determines the threshold law for re- 
combination and is, for two identical particles in general, 
an irrational number that depends only on the masses. 
Moreover, since n m i n is always a nonzero quantity, we 
find that the recombination rate is never a constant at 
threshold in 1-D. 

For bosons we have 



1023-f I 
2|<fe3| 



A23 < f (M > m) 
hz > f (M < m) ' 



1023-f I' 

while for fermions we have 



FT' 



1*2 



i 23 < f (M > m) 
623 > f (M < m) ' 



Even Parity 
Odd Parity 

Even Parity 
Odd Parity 



Note that these results do not all immediately reduce to 
the equal mass results (i.e. 023 = tt/6) of Sec. IIII Bl since 
the symmetry at 023 has not been taken into account 
here. That symmetrization eliminates some solutions, 
finally giving complete agreement with the equal mass 
results. 



IV. RECOMBINATION FOR BOSONS 
A. Cross section and event rate constant 

In the presence of a purely short-range hypcrradial po- 
tential, the scattering solution for distinguishable parti- 
cles is of the form 



d/sl 



JikK cos (0— 4>') _|_ j 



3B ' 



R 



(40) 



where the "direction" of the incident plane wave is pa- 
rameterized by the angle 0'. For 0' = 0, the incident flux 
is in the direction of the first Jacobi vector X12, while for 
0' = 7r/2, the incident flux is in the direction of X12.3. 
The quantity f 3B represents the elastic three-body scat- 
tering amplitude and is given by 



3B 



1 

2~7rfc 



E 



e im(<f>-(f>') /g2ii5„ 



1) 



(41) 



m— — 00 



We have defined the scattering phase shift in the m- 
partial wave as 5 m . The total integrated cross section 
is then [H 



Pdist 



d4\f 



3B|2 



^ m— 00 
h E 



„2i8„ 



(42) 



This expression must be modified to account for three 
separate issues. First, we are interested in the cross sec- 
tion for an inelastic process. Second, we must account for 
the appropriate identical particle symmetry, and third, 
that our asymptotic harmonics are not two-dimensional 
partial waves. 

Taking these issues into consideration, we find — fol- 
lowing the arguments of [30l ] — that the recombination 
cross section for identical bosons in terms of the S- matrix 



is 



E 



(43) 



6 



where ^^Bp denotes the final symmetrized two-body 
channel function with overall three-body parity p, and 
k labels the three-body entrance channel. At ultra- 
cold temperatures, we need only include the smallest 
n = Kmin, giving 



6 



*2£|S|k = 3 



(44) 



Similar expressions can be derived for three identical 
fermions and for systems with only two identical par- 
ticles. In fact, the only change for these other symme- 
tries — besides having different dominant k and parity 
at threshold — is the numerical prcfactor. 

We are not only interested in the cross section, but also 
in the more experimentally relevant event rate constant, 



hk 



6/v 



A", = —a= — 



/' 



^2B.3B 



(45) 



Here, we explicitly show all factors of h in order to em- 
phasize that this quantity has the appropriate units of 
length 2 /time. A3 represents the probability of a recom- 
bination event per atomic triad per unit density of atomic 
triads per unit time. The volume is the full volume 
spanned by the internal coordinates of the three parti- 
cle system, so that the density is the number of triads 
per unit two-dimensional volume. As defined in Eq. (|4"5|) . 
the rate is a function of energy. To compare with exper- 
iment, though, it should be thermally averaged to give 
the rate as a function of temperature. This averaging is 
especially important in 1-D since A3 (A) is not constant 
at threshold. 



B. Zero vs. finite-range interactions 

For identical bosons with <5-function interactions, the 
clastic atom-dimer S- matrix element is [T], 0, [2(| : 



Su = exp (2iS) = 1— -g 



24(fc 12 , 3 a) 



«(fci2,3a) 2 + 12(fci2, 3 a) - 6i 



(46) 



where £12,3 = \J^^i2.3,{E + B2). Simple analysis of 
Eq. (|46|) shows that |5ii| 2 = 1 for all energies, meaning 
that scattering in the two-body channel is always elastic. 
Hence, the amplitude for breakup (and also recombina- 
tion) is identically zero at all energies. 

For S- function interactions, the adiabatic hypersphcr- 
ical potential curves have previously been calculate d by 
solving the transcendental equation Eq. ([22]) p5l.[26l.[3l|. 
Since we want to calculate inelastic transition rates, we 
also require the nonadiabatic couplings P« and Qij in 



Eqs. (|C7[) and (|C8[) . In general, the preferred way to 
calculate these couplings numerically is from difference 
formulas. It is difficult to obtain high accuracy by differ- 
encing, though. In cases where the adiabatic solution <!> 
can be written analytic ally , however, it is possible to cal- 
culate 9$ / dR directly [32| . This is accomplished by dif- 
ferentiating the transcendental equations for k in Scc. lIIII 




FIG. 1: The adiabatic potential curves (a) and first deriva- 
tive couplings (b) for 5-function interactions (a = 2). The 
adiabatic potential curves (c) and first derivative couplings 
(d) for the Poschl- Teller potential with L = a = 2. All axes 
are dimensionless. 



to determine equations for On/ dR. Rather lengthy ex- 
pressions for Pij and Qy result, but they do allow the 
couplings to be calculated essentially exactly. The first 
few elements of the first row of the antisymmetric matrix 
P are shown in Fig. [IJb) and (d). The important fea- 
ture of this figure is simply that these elements are not 
zero. They couple the two and three-body channels, yet 
we know from Eq. (|46[) that the amplitude for any pro- 
cess connecting the two and three-body channels must 
vanish. We will numerically demonstrate that despite 
the nonzero couplings, the solution to the coupled radial 
equations indeed gives vanishing probability for inelastic 
events. 

In order to facilitate this demonstration, we consider 
the Poschl- Teller two-body potential, 



V(x) 



-D sech 2 — , 



D > 0, 



(47) 



which give s a Schrodinger equation having an analytic so- 
lution [331 ]. Defining 7/ = vl-j- ADmL 2 , the zero-energy 
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scattering length for the even parity solution is 



2 V 



tt sec rj— 



(48) 



In this expression, H x = Jem + 4>(x + !)■ 7EM is the 
Euler-Maschcroni constant, and ip(x) = r'(x) /T(x) is the 
digamma function. The energy eigenvalues are 



E 



(2n + 1- 77) 
4mX 2 



0,1,2. 



(49) 



The scattering length becomes infinite when an even par- 
ity state sits at threshold, that is, when E = for 
n = 0,2,4.... 

To check whether recombination with the Poschl- Teller 
potential recovers the S- function result (|46|) , we set m = 
1, fix the scattering length to a = 2, and consider the 
effect of letting L — ► 0, approaching the 5-function limit. 



When L = 2 and D 



2 ■ 



the Poschl- Teller potential 



gives the same two-body binding energy and scattering 
length as the 5-function. The only difference is that the 
former is of finite range. The difference in the resulting 
three-body potential curves and couplings, however, is 
more subtle. These potentials and couplings arc shown 
in Fig. [2(c) and (d). By comparison with the S- function 
results in Fig. [Ha) an d (b) , we see that all features for 
the Poschl- Teller results tend to be pushed to larger R. 

Our numerically obtained |Sii| 2 for both kinds of po- 
tentials are shown in Fig. [2l Below the dimer break- 
up threshold at E = 0, the collision is purely elastic. 
Above the break-up threshold, the phase-shift in the elas- 
tic channel acquires an imaginary part, which appears as 
a deviation from |<Sii| 2 = 1, implying a nonzero probabil- 
ity for break-up of the dimer. As indicated in the figure, 
our numerical calculations show that as L — > 0, the colli- 
sions do indeed become purely elastic at all energies, in 
agreement with Eq. (|46[) . Since the S'-matrix is symmet- 
ric under time-reversal, the amplitude for recombination 
also vanishes. It should be stressed that the solid black 
line in Fig. [5] is the result of our numerical calculation 
for the (5-function, and not simply a plot the analytical 
result, Eq. (|46|) . 
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FIG. 2: The convergence of ISnl 2 for finite -range potentials 
(with range L) to the zero-range theory is shown. The scat- 
tering length is held fixed at a = 2 as L in Eq. (|47[) is taken to 
zero. The data for L = 1/8 (squares) are nearly indistinguish- 
able from the solid 5-function line. All axes are dimensionless. 



the limit a 
take 



00. We follow the insight of Rcf. [34| and 



V, 



A 2 



eff\ 



[c- 2d + 4:d(Ax) 2 ] cxp(-A 2 a- 2 ). (50) 



The prefactor A 2 /m gives the interaction units of energy 
in atomic units, and m is the mass of the helium atom 
in atomic units (m = 7296.299m e ). One can tune the 
free couplings c, and d so that the potential reproduces 
physically reasonable values for the scattering length and 



effective range of, for example, He [3f 



208 a.u. and 



ro = 14 a.u.. A reasonable requirement of any renormal- 
ized model is that low-energy observables remain inde- 
pendent of the momentum cutoff A. We find that the 
values c = 1.97151 and d = -1.46165 for A = 0.16 a.u. 
yield the desired effective range parameters. These values 
also give a two-body binding B2 = 3.39985 x 10~ 9 a.u.. 
We have verified that B 2 varies by less than one percent 
over a wide range of cutoffs 0.08 < A < 0.3, and that the 
couplings c and d are of order unity over this range. 



D. Low-energy scaling behavior 



C. Low-energy effective interaction 

While the potential Eq. (|47|) has the advantage of yield- 
ing an exact analytic two-body solution, it does not facil- 
itate the study of recombination into a single high-lying 
two-body state in the limit a — > 00. The disadvantage 
stems from the fact that the potential is purely attrac- 
tive. As is well known, any purely attractive potential in 
1-D will support an even parity bound state no matter 
how small the coupling. Hence, for the potential Eq. (|4"T]) 
to have a — > 00 with a single bound state, we must let 
D — > 0. Therefore, we are motivated to construct a 
renormalized low-energy potential model that supports 
a single (shallow) bound state with finite couplings in 



From Eq. (fT4")) and Sec. IIII[ we know that the inelastic 
transition probability must scale as IS12I 2 ~ k e ~ E 3 
at threshold 3 for three bosons. This be- 

havior is already plausibly seen in Fig. [2] since this scal- 
ing for S12 implies 1 — |<Sii| 2 ~ E 3 at the three-body 
breakup threshold. When the scattering length a is the 
longest length scale, we also know from Eq. (|15p that 
1 5*12 1 2 ~ (ka) 6 ~ E 3 a 6 . We demonstrate this behavior 
quantitatively in Fig. [3J The points in the figure were 
obtained by first tuning the effective potential Eq. (|5"0|) 
to reproduce a given scattering length, and then cal- 
culating the recombination probability near threshold, 
E — > 0. Finally, we divided the probability by E 3 to ex- 
tract the constant of proportionality. The constant A$ 
in I S12 1 2 = AqE 3 is plotted as a function of the scatter- 
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FIG. 3: A numerical demonstration of a scaling of the three- 
body recombination probability is shown. The black line is 
Aq = 3.0 x 10 9 a 6 , and the circles are our numerical results. 
The units of Ao are atomic units of inverse energy cubed con- 
sistent with the use of Eq. [50] to model the 4 He two-body 
interaction. 



with 



12 



7o 



3 1 /42V2 7r 



-2.052277. 



(53) 



We thus find an attractive, universal potential when 
a > 0, and a similarly universal, repulsive potential 
when a < 0. In the former case, this universal potential 
converges to the highest-lying two-body threshold while 
for the latter, it is the lowest potential converging to the 
three-body breakup threshold. 

This result can be derived in several ways. For in- 
stance, the (5-function transcendental equations Eqs. (fH?|) 
and (f2"2"|) for positive and negative scattering lengths, re- 
spectively, can be solved in the small q or n limit. The ex- 
act same universal curves can similarly be extracted from 
the analogous quantization equations for the square- well 
interaction [Eq. ([2TJ)) . for example] using 



ing length in Fig. [31 The a 6 scaling is clearly seen when 
compared to the solid line which is Ao = 3.0 x 10 9 a 6 . 

It is interesting to note that the lack of inelastic pro- 
cesses for zero-range interactions is actually a conse- 
quence of perfect destructive interference in the exit 
channel. Indeed, since the couplings Pjj and Qij are non- 
zero, the only way for the inelastic probability to vanish 
is through some sort of interference effect. It is possi- 
ble to demonstrate this perfect interference by adding 
an arbitrary short-range three-body interaction Vs(R) of 
characteristic length L% to the zero-range hyperradial po- 
tentials. The short-range three-body interaction destroys 
the perfect interference and leads to a nonzero recombi- 

1— is l 2 

nation probability. Considering the ratio C — / ka \l 
near threshold, we find a nonuniversal power law (that 
depends on the short-range nature of V3) of the form 
C ~ (L^/aY 1 . Although the particular value of c\ is 
nonuniversal, we always find c% > such that K 3 — > as 



V. UNIVERSALITY IN ONE DIMENSION 

In three dimensions, the three-body problem exhibits 
universal features in the limit \a\ —> 00 related to Efimov 
physics [3, [l6|, [3(|. Here, we consider this limit in 
one dimension. For L <C R <C \a\ — with L generally the 
characteristic size of the two-body interaction — we find 
that 



YluR 



(51) 



to leading order in R/a for a > and a < 0, respectively 
The corresponding adiabatic potentials are 



U(R) 



7o 
2[iaR 



(52) 



L 



n > 0, 



(54) 



assuming there are n two-boson bound states. This 
choice for (3 follows from the fact that the two- 
body square-well has a zero-energy bound state when 
\Z2Jp2bVqL = mr. Reducing this phase slightly gives 
a < 0; and increasing it, a > 0. Finally, the universal 
curves can be obtained quite generally from the quan- 
tization conditions for arbitrary short-range potentials, 
Eq. (|33[) . using this expression for the log-derivative in 
the limit 2^2bEL 2 < 1: 



1 



(55) 



It is clear that the universality of the 1-D problem has a 
different character than the 3-D problem. First, 70 is the 
coefficient of R^ 1 instead of R~ 2 as in 3-D. Second, the 
universal potential Eq. ([B"2"l) actually vanishes in the limit 
\a\ — > 00. The lowest three-body continuum potential 
in this case is instead an attractive R~ 3 potential. For 
square-well interactions, we can derive it explicitly using 
(3 from Eq. JM]) with \a\ — ► 00 in Eq. ([%?])■ yielding 



(56) 



This potential, however, is not universal, although nu- 
merical calculations with four different short-range two- 
body potentials suggest that two aspects are general: (i) 
the attractive R~ 3 behavior and (ii) the increasing inter- 
action strength with n. 

As \a\ — > 00, the solution in the lowest channel in the 
region L -C R <S \a\ will obey the equation 



1 d 



R 







2nRdR\ dRJ 2^aR 



la 



Qu{R) 



F(R) = EF(R). 

(57) 

The diagonal nonadiabatic coupling —Q\\{R) is always 
repulsive and falls off as Rr 2 , while the universal 
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FIG. 4: The lowest adiabatic hyperspherical potential curve 
for the two-body interaction in Eq. (|50[) with several different 
values of a < (a given in a.u.). The curve is multiplied by 
the factor 2fiRa to give a unitless quantity and more clearly 
reveal where the potential behaves as a repulsive 7i _1 . The 
horizontal solid line is the universal number 70 ~ —2.05 (|53l) . 
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term will vary smoothly from an attractive to a repulsive 
potential as a varies from +00 to —00. When \a\ = 00, of 
course, the i?" 1 potential vanishes, leaving just the R~ 3 
potential in Eq. (|5fJ|) . 

In order to demonstrate that 70 is indeed a univer- 
sal constant, we again turn to the effective potential 
Eq. . We have calculated three-body potential curves 
Eq. (|50p giving different two-body scattering lengths (the 
effective range is held constant at 14 a.u.). Figurc[4]shows 
the lowest three-body potential for increasingly negative 
scattering lengths. The potential curves have been mul- 
tiplied by the factor 2(iaR in order to more clearly reveal 
the universal behavior. We see that the curves do in fact 
approach 70 over a range of R consistent with the con- 
dition ro^Ci?<C|a|. Ati?^>|a|, the potentials again 
approach the three-body breakup threshold with the R~ 2 
behavior predicted in Sec. IIIII (which translates to 
as plotted in the figure). 

For large positive scattering lengths, the lowest poten- 
tial curve — the two-body channel — supports a three- 
body bound state so long as there is at most one weakly 
bound two-body state. Figure [3] shows the lowest poten- 
tial curve for various values of the two-body scattering 
length along with the hyperradial wave functions for the 
three-body bound state in this channel. It is clear that 
the universal R~ x portion of the curve has a strong in- 
fluence on these states since a significant portion of the 
necessary phase is accumulated in this region. This ar- 
gument is supported by a simple WKB calculation using 
the the Bohr-Sommerfeld quantization condition 

S2 dR h E -^ = { n+1 ^ (58) 

for classical turning points Ri. We note that the Langer 
correction does not appear in this equation since it can- 
cels the attractive R~ 2 term obtained from eliminating 
the first derivative in the hyperradial kinetic energy — 
which is required to use this WKB phase integral. The 
energies of the nodeless solutions as a — > 00 using the 
above equation are in reasonable agreement with a B- 



FIG. 5: (a) The lowest adiabatic hyperspherical potential 
curve for the two-body interaction in Eq. (|50[) with several 
different values of a > (units of a are a.u.). The curve is 
multiplied by the factor 2(iRa to reveal the region in which 
the potential behaves as an attractive J? -1 . The horizontal 
solid line is the universal number 70 ~ —2.05 (|53[) . (b) The 
hyperradial wave functions for each potential in (a). The wave 
functions have been scaled such that their maximum is unity. 
As a — > 00, we expect the three-body binding B3 — > 4E?2, as 
is the case for 5-function interactions 



TABLE I: The three-body bound-state energies are shown 
in atomic units for various two-body potentials of the form 
Eq. (|50[) tuned to give the scattering lengths in the first col- 
umn. The WKB estimates improve as a — > 00. 



a 


Numerical B-spline 


WKB 




50 


-1.8644 x 10~ 7 


-4.8391 x 10" 


-8 


800 


-6.0569 x lO" 10 


-4.8839 x 10" 


1(1 


12800 


-2.5156 x 10~ 12 


-2.6829 x 10" 


12 


204800 


-1.0324 x 10~ 14 


-1.1533 x 10" 


1-1 



spline calculation using the numerical potential curves. 
These results are tabulated in Table HI 

Further examination of Eq. (|58|) also shows why — de- 
spite the long-range R~ 1 behavior — there is only a single 
three-body bound state. If we evaluate the WKB phase 
at the two-body threshold for R between L and a, then 
we find that it lies between ir/2 and 0.8177T for all a/L 
between 8.07286 and infinity. The phase remains finite 
because the strength of the R^ 1 potential decreases with 
a at the same time that the domain over which it holds 
grows with a. So, there is sufficient phase accumulated 
in this universal region alone to support a single bound 
state by Eq. (|58|) . Moreover, the phase contributed from 
the small- R region, R < L, would have to exceed roughly 
0.77T to produce a second bound state. While this is not 
impossible, it does not seem likely. 
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turns out that odd parity, identical fcrmions share this 
threshold law in the limit that the two-body odd-parity 
partial wave scattering length — appropriate for fermion- 
fermion interactions — goes to infinity. The threshold 
scaling for odd parity bosons and even parity fermions 
is also changed from the predictions of Sec. In the 
limit that their respective scattering lengths arc infinite. 
K 3 oc E 3 for both cases. 



FIG. 6: The inverse of the atom-dimer scattering length mul- 
tiplied by the atom-atom scattering length is shown as a func- 
tion of the inverse of the atom-atom scattering length. In the 
limit that a — > oo, we find 



Finally, we consider the atom-dimer scattering length 
a ad in the universal limit. Again using Eq. l50l we calcu- 
late a ad as a — > oo. Our results (plotted in Fig. [6]) show 
a clear trend towards clad ~ 10a. This finding is con- 
sistent with the presence of a bound-state at the atom- 
dimer threshold in the 5-function model [l], HH, [3l[ ■ As 
a — > oo this zero-energy resonance results in a ad — > oo. 

Before we end our discussion of universality, it is worth 
mentioning the consequences for the recombination rate. 
So long as a is finite, the lowest three-body continuum 
potential still behaves as predicted in Sec.|TTT]for R^> \a\. 
The power-law scaling of the S'-matrix with E and a is 
thus the same as found in that section. The a-dependence 
is modified, though, by the universal region of the poten- 
tial. Using arguments similar to those in Refs. [l7l. [37|. 
we can use WKB to determine the modifications. For 
a < 0, the coupling driving recombination peaks around 
R = L, so the system must tunnel through the repul- 
sive R~ 2 potential at R ^> \a\ and through the repulsive 
universal in the region L <C R <C |a| in order to 
recombine. The WKB tunneling integral then leads to 
the modified threshold scaling for a < 



A3 oc — (fca) 6 exp 




(59) 



where £ is a numerical constant on the order of unity that 
depends on the exact range of R over which the universal 
potential is valid. 

For a > 0, in analogy to the 3-D case, recombination 
can be modified by interference to give 



K 3 oc -(fca) 6 sin 2 




(60) 



where $ is the short-range, fid, phase accumulated 
in the non-universal portion of the two-body channel. 

When \a\ = 00 and the incident three-body continuum 
potential is given by Eq. (f56|) . the threshold scaling is 
nontrivially modified. The attractive R~ 3 potential is 
equivalent to K min = 0, so the recombination rate is, in 
fact, independent of energy at threshold, A3 oc const. It 



VI. DISCUSSION AND SUMMARY 

The work presented here has been carried out in 
strictly one dimension. It is worth commenting on the 
relation of this study to experimentally-realizable, effec- 
tive 1-D geometries. Olshanii [lj| has determined the ef- 
fective 1-D scattering amplitude for two particles with 3- 
D s-wave scattering length a^o interacting under strong 
cylindrical harmonic confinement. His analysis leads to 
the following one-dimensional effective range expansion: 



fc tan 5 = — 
a 



C(3/2)ai 1 



8a 2 



2 



+ 0(fc 4 ), (61) 



where a± is the oscillator length in the confined direc- 
tion, and is the Riemann zeta function. The 1-D 
scattering length is thus determined at zeroth order in k 
to be 



2a 3 _o 



l + C(l/2) 



0-3D 



(62) 



One may argue that the 3-D effective range should also 
be present in the fc 2 term of Eq. (foTj) , but since a± ^> Tq D 
in general, it is reasonable to assume that the contribu- 
tion from Tq D is small compared to a±. With this ap- 
proximation in hand, it is possible to calculate the 3-D 
parameters a± and a^u that correspond to the 1-D pa- 
rameters a and tq that we have used. For r$ = 14 a.u. 
and our larger values of a, we find experimentally acces- 
sible values for a± and ol^d- All values are tabulated in 
Table [TTJ One notable point is that the 3-D scattering 
length is negative for positive 1-D scattering lengths. So 
while there is no shallow two-body bound-state in the 
3-D system, a shallow bound-state appears as a result of 
the cylindrical confinement. 

Finally, we comment on the role of three-body recom- 
bination in one recent experiment. The experiment by 
Tolra et al [1] has measured three-body recombination 
rates in order to probe properties of the many-body wave- 
function. This idea was originally proposed by Kagan et 
al 38], who showed that the event rate constant K3 is 
proportional to the three-body local correlation function 
#3. Therefore a measured reduction in A'3 from a 3-D 
system to a 1-D system in Rcf. Q was interpreted as a 
reduction in g% and a clear signature of enhanced corre- 
lations. 

While a complete description of three-body recombina- 
tion in atomic waveguides requires a full-scale calculation 
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TABLE II: The 3-D parameters a^o and a± are calculated to 
give the desired 1-D effective range rg — 14 a.u. and the 1-D 
scattering length in the first column. 



a (a.u.) 


a-iD (a.u.) 


a± (a.u.) 


50 


-73.66 


47.50 


200 


-63.62 


119.70 


800 


-78.46 


301.62 


12800 


-160.84 


1915.14 


204800 


-377.39 


12160.40 


3276800 


-925.65 


77213.80 



of three particles in a confinement potential, we attempt 
to address some of the issues within our 1-D framework. 
Our work suggests an alternative explanation of the ob- 
served suppression of K3 in terms of the three-body hy- 
perradial potentials. For the experimental parameters of 
Ref. 3 (a 3D ~ 100 a.u., a± w 1100 a.u.), Olshanii's for- 
mula Eq. (frJTj) gives a large negative 1-D scattering length 
a w —5300 a.u. indicating the absence of a shallow 1-D 
bound state. The wavefunction in the lowest three-body 
entrance channel must therefore tunnel under both the 
repulsive n^ nin /{2^R 2 ) potential in the region R ^> \a\ 
and the repulsive 7o/(2//ai2) (recall that 70 and a are 
both negative) potential in the region a± <C R <C |o| 
before reaching the region R < a± where a recombi- 
nation event into a deep two-body channel may occur. 
Therefore the measured suppression in the confined ge- 
ometry could be the combined effect of the threshold law 
K\ D oc ^-(fca) 6 and the universal barrier given in 
Eq. l[5"2"]) leading to the suppression given in Eq. 

While we are confident that we have solved the 1- 
D Schrodinger equation accurately, it is difficult to 
make quantitative claims regarding the measured sup- 
pression Q. One difficulty arises from the fact that the 
temperature of the system in Ref. @ is not very well 
characterized, and therefore, it is unclear what energy 
three-body collisions occur at. In addition, there are a 
number of theoretical complications which can only be 
accounted for by performing a full-scale calculation of 
three interacting bosons in a confined geometry. For in- 
stance, since there are excited radial modes of the trap 
and the two-body binding energy is typically much larger 
than the spacing of these radial modes, there will be a 
series of two-body thresholds attached to each excited ra- 
dial mode. All of these thresholds arc open at the initial 
three-body threshold energy, and our calculations have 
not accounted for this. Also, our strictly 1-D calcula- 
tions undoubtedly do not correctly represent the hyper- 
radial potentials in the region R < a± where the system 
is neither purely 3-D nor strictly 1-D in nature. Since the 
potential appears in the exponent of the WKB tunneling 
integral, any uncertainty in the potential leads quanti- 
tatively to very different results for the suppression of 
K%. Nevertheless, in view of our findings regarding uni- 



versality in 1-D systems, it is possible that the observed 
suppression of K 3 in Ref. [|| was a direct probe of the 
universal potential Eq. l[5"2"j). 

We have shown in this study how inelastic processes 
such as three-body recombination and collision induced 
break-up behave near threshold for all combinations of 
identical particles. For the system of three identical 
bosons, we have further investigated the behavior of these 
processes with respect to the range of the pair-wise in- 
teractions and to the two-body scattering length. Wc 
have demonstrated numerically that the probability for 
inelastic events vanishes as the range of the interaction 
is taken to zero in agreement with previous analytic re- 
sults [l|. Our work is related to a recent Letter that 
has shown how inelastic processes could in fact be possi- 
ble if one considers a zero-range two-channel model with 
a confinement induced Feshbach resonance [3^]. Finally 
we have explored in detail the nature of universality for 
three-body systems in 1-D. 
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APPENDIX A: SYMMETRY OPERATIONS 

The pair permutation operators P,j have the following 
effects on the hyperangle: 

Pl20 = 7T - <j) 
7T 

P 3 l<t> = y - 4>- ( A1 ) 

Solutions of definite parity may be found by via the op- 
eration 

LL/> = 7T + 4>. (A2) 

For three identical bosons, the symmetric projection re- 
quires 

S = 1 + P12 + P23 + P31 + P12P23 + P12P31 ■ (A3) 

Similarly, for three identical fermions, the antisymmetric 
projection operator is 

A=l-P 12 - P23 - P31 + P12P23 + P12P31 ■ (A4) 

To determine the boundary conditions required by per- 
mutation and parity, it is useful to apply these operators 
to the free particle solutions of Eq. (fTI))) (l8j . Doing so, 
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we find that the simultaneous symmetry eigcnstates are 
real harmonics with indices that are multiples of three. 
Explicitly, for bosons 



cos 3ncA, 
sin 3n0, 



and for fcrmions, 

sin 3n</>, 



oc 



cos Zn<j) 1 



n 

n - 



n - 
n 



0,2,4,. 
1,3,5,. 



2,4,6,. 
1,3,5,. 



Even Parity 
Odd Parity 



(A5) 



Even Parity 
Odd Parity 



(A6) 

From these expression, we conclude first of all that we 
can reduce the integration range by a factor of 12 - 
as expected — from 27r to only 7r/6 in cf>. We can also 
conclude, for instance, that three identical fermions in 
the odd parity state must obey the boundary conditions 
$'(0) = and $(tt/6) = 0. 

If there are only two identical particles, then we need 
the operators 5 = 1 + P12 for identical bosons and 
A = 1 — P12 for identical fermions. The simultaneous 
symmetry eigcnstates for bosons are now 



cos ntp, 
sin no 




n : 



72 = 

n 



= 0,2,4,. 
1,3,5,.. 

2,4,6,.. 
: 1,3,5,. 



Even Parity 
Odd Parity 

Even Parity 
Odd Parity 



The range of integration can in this case be reduced by 
a factor of 4, and the boundary conditions extracted as 
above. 



We begin with the variational expression 



E 



J d<j)RdR $(R, <f>)H(R, cA)*(i?, , 
Jdcj)RdR |^(i?,cA)| 2 



(CI) 



assuming "J is real without loss of generality. The Hamil- 
tonian is given in Eqs. ([8]) and (J9j> . Integrating the hy- 
perradial kinetic energy by parts in Eq. (|Clj) gives the 
variational expression, 



b = 



Jg d4>RdR [-|*|* + Vjk* - 2^)*) 
J d<t) [R 2 \^(R 2 ,^+R 1 \^(R 1 , ( f>)\ 2 } 



(C2) 



where b is minus the log-derivative [4Ct |44| normal to a 
hypersphere, 



din* 

dh 



(C3) 



APPENDIX B: SUMMARY OF RESULTS FOR 
SQUARE- WELL AND ^-FUNCTION 
INTERACTIONS 

We summarize in this section the boundary and quan- 
tization conditions, along with the large R solutions, for 
all symmetries of three identical particles with square- 
well interactions (see Table IIIII and Sec. MI Bp and for 
all symmetries of two identical particles with (^-function 
interactions (see Table ITVl and Sec. IIII Cj) . 

We expect the square- well results in Table IIIII can be 
generalized to an arbitrary short-range interactions in ba- 
sically the same way as Eq. ([33]) . That is, the left-hand 
sides of the quantization conditions would get replaced 
by the two-body log-derivative. 



In the adiabatic hypersphcrical representation, 



9{R, <t>) = Fi(R)$i(R; 4>) (C4) 



where $i(i?;</>) are eigenstates of Eq. © with eigenval- 
ues Ui(R). Substitution of the expansion Eq. (|C4[) into 
Eq. (|C2[) results in a set of coupled equations for Fi(R). 
We further expand the radial functions in a basis set as 
Fi(R) = J2 n c mB n (R), converting the coupled equations 
into the generalized eigenvalue problem 



APPENDIX C: R-MATRIX PROPAGATION 

In order to solve for the multichannel S'-matrix, we use 
the eigenchannel i?-matrix method coupled with the adi- 
abatic hyperspherical representation. The development 
closely follows Refs. and [iH, S3, S3] - 

TABLE III: The boundary conditions, quantization condi- 
tions, and allowed asymptotic k for three identical particles 
with square-well interactions. The constants ft and (f>b are 
defined in Sec. HITBl 



bAc = Tc. (C5) 



The matrices A and T arc defined as 



Even Parity Bosons 
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TABLE IV: The boundary conditions, quantization conditions and large R solutions are shown for the two-identical-particle 
case. The constants as, Qr>, and (j>23 are defined below Eq. (|34[1 and in Appendix HT1 



Even Parity Bosons 




$'(0) = and $'(f + e) = -$'(§ - e) 




(— ft + ly 2 g 2 otsot. D R?) tan [ft(0 2 3 - f )] - /jgas-Rtan k0 2 3 tan [ft(</> 2 3 - §)] + fttanK0 2 3 = 


Atfl(as + 2a D )R 


K = rm/\(f>23 - ||, n = l,2,3... and ft = (n + 1/2)tt/|<^ 2 3 | , n = 0, 1, 2. 




Odd Parity Bosons 




$(0)=0 and $'(f +e) = -$'(§ - e) 




(— ft + -|^i 2 (; 2 asQ:D-R 2 ) tan [k(</!>23 — § )] — figasR cot K0 2 3 tan [ft(023 — f )] + ft cot K023 = 


fig(as + 2oid)R 


K = nir/\(f>23 — § | 5 ti = 1,2, 3... and ft = n7r/|023], n = 0,1,2... 




Even Parity Fermions 




$(0) = and $(f ) = 




ft COt [ft(023 — f )] — ftCOtft023 = l^gUoR 




ft = n7r/|<^23 — f |, ?i=l,2,3... and ft = n7r/|0 2 3|, n = l,2, 3... 




Odd Parity Fermions 




$'(0)=0 and $(§)=0 




ftcot [ft(</>23 — f )] + fttanft023 = Ip^ganR 




ft = rar/|023 — f |, n = 1,2, 3... and ft = (n + |)7r/ 023 1, « = 0, 1, 2.. 





A(m,i),(nj) — ^i.j/J [RlB n (Ri)B m (R\) + R2B n (R2)B m (R2)} 



(C6) 



(m,i),(raj) 



dB m 8B n n /" <9£? m 95 



Hi 



V%B m B n + S itj B m [k 2 - 2fxUi(R)] B n 



We require the first-derivative nonadiabatic coupling 

d 



P ij (R) = (( MR 



dR 



and the second derivative coupling, 

d 2 



Qij(R) 



MR; 



dR 2 



(C7) 



(C8) 



The double bracket notation implies integration only over 
(/). The square of P is given by the symmetric part of Q 

dPjj . 
dR : 



through the relation Qij 

'OMR 



p2 



dR 



d*j(R);9 



dR 



(C9) 



Equation (|C5|) is solved for the expansion coefficients 
c, yielding a solution in the region R\ < R < R2. For the 
first region < R < Ri, we require F((0) = but impose 
no boundary condition at R\ . The generalized eigenvalue 
problem then yields an eigenvalue bp and wave function 
Fip for each open or weakly closed channel at R\ . 

To propagate the solution to large R, we solve the sys- 
tem of equations again in the region Ri < R < R 2 with 
no boundary conditions at cither R\ or i?2- 



A valid solution in the region < R < R2 is then con- 
structed by matching the solutions from the two regions 
at R\ and by requiring that the overall solution be an 
eigenchannel solution of the i?-matrix, 



R= (F)(F')~ 1 . 



(CIO) 



That is, ^ should have constant normal derivative at the 
surface R = R2 . This procedure is repeated until the so- 
lutions can be accurately matched to analytic asymptotic 
forms HUES!- 

To be more explicit, we write the full wave function 
outside the i?-matrix volume as 



V°a ut (R,< 



3 



MR; 0) [f^R^-g^R)^] , (Cll) 
where fj(R) and gj(R) are 



, 2 „ cos krR, if 7 is a two-body channel 



J K (kR), 



if j is a three-body channel 



/ 2 5 sinkfR, if 7 is a two-body channel 
gj(R) = {V* k t a 5 

\Y K (kR\ if j is a three-body channel 

(C12) 
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The order n is determined as described in Sec. IIIII The 
solution inside the volume involves the numerical func- 
tions Fjp 



¥?(R;<t>) = Y j F j 0(R)$ j (R; 



3 



(C13) 



The matrices I and J are obtained from matching at 
some large distance R = Rm, which is conveniently ac- 
complished using 



i jP = w( Sj ,F jfi yw( Sj ,f j ) 



(C14) 



where W(f, g) denotes the Wronskian of / and g. Defin- 
ing (f). ti = Sijfi and (g) y = Sijgi, the if-matrix is 



K(i? M ) = (f-f'R)(g-g'R)- 



(C15) 



The i?-matrix is calculated via Eq. (|C10[) using the 
numerical solutions F^ and the result F^(Rm) = 
—bpFip(RM)- The latter holds only at large R where 



Pij — > since the exact relation is 



bpFi^Flp+YjPijF, 



(C16) 



in which each quantity is evaluated at Rm ■ This relation 
is obtained by differentiating Eq. (|C13[) and projecting 
the result onto <5j. Finally, the S- matrix is calculated 
from K using 



S(Rm) = 



t + iK(R 



Mi 



1 - iK(R 



M 



(C17) 



It is important to propagate the i?-matrix to large R in 
order to obtain a converged, unitary S'-matrix. In Fig. [71 
we show the convergence of a few S'-matrix elements with 
the matching distance Rm for an eight-channel calcula- 
tion using the Pdschl- Teller potential Eq. (|4"Tj) with L = 2 
and D = 1/2. Figure [7fa) shows the lowest nine poten- 
tial curves corresponding to three-body channels. (Since 
it converges to -1 on the scale of the figure, the two- 
body channel is not visible.) The horizontal dashed line 
shows the energy at which the calculation in Fig. |7fb) 
was carried out, and the vertical dotted lines mark the 
classical turning points for the first three channels. Note 



that the probability \Su(Rm)\ peaks approximately at 
the classical turning point for Ui(R). 

For the calculations presented in Fig. we propagated 
the ii-matrix to R — 1000a to assure convergence. It 
is evident from Fig. [7] that the quantity IS12I 2 has al- 
ready converged by R — 200a. We have also verified that 
this probability is stable with respect to the inclusion of 
more coupled channels. The other probabilities, however, 
have not yet converged as well, although their magnitude 
makes them negligible. 

Our calculation of the matrix elements in Eq. (|C6|) 
is simplified by using B-splines as our radial basis set 

°- 8 l I: ' :\l : ' IV \ \ I 1 I 1 I 1 I 1 I 1 I 1 T 

(a) 




FIG. 7: (color online) (a) Three-body potential curves for 
the Poschl- Teller two-body potential with L = 2 and D — 
1/2. The horizontal dashed line indicates the collision energy, 
and the vertical dotted lines mark the classical turning point, 
(b) The convergence of the transition probability between the 
two-body channel and the lowest three three-body channels 
at E = 0.4£>2 as a function of the matching distance. 



{B n (R}} [43. This choice also simplifies the imposition 
of boundary conditions since B-splines have only local 
support. We typically use ten fifth-order B-splincs within 
each i?-matrix sector, leading to a (10 x N) x (10 x N) 
matrix equation (N is the number of channels). The size 
of the sectors R4+1 — Ri is chosen to be no more than one 
de Broglie wavelength in the lowest (two-body) channel. 
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